Finite-time scaling of dynamic quantum criticality 
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Time plays a fundamental role in quantum criticality owing to the interplay of static and dynamic 
behaviors. Yet, it is a great challenge to study and understand real-time quantum critical dynamics 
both at zero temperature and at finite temperatures. Here, we consider the competition among an 
external time scale, an intrinsic reaction time scale and an imaginary time scale arising respectively 
from an external driving field, the fluctuations of the competing orders and thermal fluctuations and 
explore a regime of finite-time scaling. Through a successful application in the dark impulse regime 
of quantum Kibble- Zurek mechanism at zero temperature and the first solution of real-time Lindblad 
master equation near a quantum critical point at nonzero temperatures, we show that finite-time 
scaling offers not only an amenable and systematic approach to detect and determine dynamic critical 
properties, but also a unified framework to understand and explore the nonequilibrium dynamics of 
quantum criticality. 

PACS numbers: 64.70.Tg, 64.60.Ht, 75.10.Pq 



Detecting and understanding dynamics of quantum 
phase transitions (QPTs) are of great challenge and im- 
portance [Il-IH . Recent experimental breakthrough in 
ultracold atoms @ promises new tools to study quan- 
tum critical dynamics Different from classical phase 
transitions whose static and dynamic fluctuations are 
independent, "time" plays a fundamental role in quan- 
tum criticality owing to the interplay of static and dy- 
namic behaviors. A continuous quantum phase transition 
(QPT) taking place at zero temperature is dominated by 
a reaction time r s , which arises from the fluctuations of 
the competing orders and blows up as r s ~ \g\~ vz with 
the standard critical exponents v and z as the distance to 
the critical point g vanishes [l| . It describes how fast the 
system reacts to external perturbations. A protocol 
in which the effect of r s is clearly exhibited is to drive 
a system through its critical point by increasing g, say, 
with a time rate of R. In this case, there exits a charac- 
teristic time scale associated with R . This is a driving 
time scale Td ~ R~ z / r with a rate critical exponent r 
that is related to z and other static critical exponents 
depending on which controlling variable is being swept 
|20j. Consequently, no matter how small R is, the adia- 
batic evolution for r s < breaks down at a freeze-out 
time instant — t and resumes at t with t ~ r^. In this 
adiabatic-impulse-adiabatic approximation of the quan- 
tum Kibble-Zurek mechanism (KZM) j8-ll|, the condi- 
tion for the breakdown of adiabaticity was employed to 
reckon a famous KZ scaling (see below); while the state 
at —t is assumed to be frozen until t, at which it is re- 
activated and serves as the initial state for the subsequent 
evolution. Apparently, however, this is not the real situa- 
tion, as for example, dynamical scaling has been reported 
within this impulse regime [l3, 14 1. 

Natural systems and their measurements exist in- 
evitably in nonzero temperatures, though probably only 
initial thermal states need considering in ultracold atoms 



The finite temperature T introduces an "imagi- 
nary" time scale tt = 1/T (in units in which the Plank 
and the Boltzmann constants are set to 1) since the real 
time is its analytical continuation to imaginary num- 
bers through a Feynman path integral representation [l[ . 
Whereas for tt > t s , thermal effects can be regarded as 
perturbation to the quantum ground states and conven- 
tional quasiparticle approaches are usually effective; for 
tt < T s , thermal fluctuations dominate. Yet, it is this 
latter quantum critical regime (QCR) that is suggested 
to the arena of exotic behavior exhibited in a wide range 
of strongly correlated systems [D-d, [l2[ . As both phases 
exhibit complex long-range quantum entanglement ow- 
ing to the long t s and are violently excited thermally in 
the QCR, new approaches have to be invoked. 

Here we systematically study quantum critical behav- 
ior due to the interplay among the three time scales r s , 
Td and tt according to the theory of finite-time scal- 
ing (FTS) [Hi Hj. As seen in Fig. [TJ besides the usual 
equilibrium regimes and the QCR which are respectively 
dominated by t s and tt, the most important consequence 
is that a new nonequilibrium FTS domain is created. In 
the impulse regime of the KZM at T — 0, although the 
system falls out off equilibrium, the state does not cease 
evolving; rather, it evolves according to the imposed time 
scale Td instead of r s . As a consequence, nonadiabatic ex- 
citations obey FTS just within the dark impulse regime, 
which overlaps the FTS regime. Therefore, FTS improves 
the understanding of KZM on its dark impulse regime. 
Similarly, at T =fi 0, in the FTS regime, there are now 
nonadiabatic thermal excitations controlled again by 
and thus obeying again FTS. This then provides an extra 
excess to detect and study quantum criticality at finite 
temperatures. Owing to its conceptual simplicity and its 
accessability, FTS therefore provides a unified framework 
not only to determine dynamic critical properties, but 
also to understand and explore the nonequilibrium dy- 
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namics of quantum criticality both at T = and T ^ 0. 

In the following, we shall first develop the theory of 
FTS for QPTs. Then we shall exemplify the determi- 
nation of the critical point and critical exponents in the 
impulse regime at T — using the paradigmatic quantum 
Ising model. Finally we shall verify FTS by solving for 
the first time the Lindblad master equation 2l], 22 1 near 
the quantum critical point at finite temperatures, which 
also shows that this equation is an appropriate platform 
for studying real-time dynamic quantum criticality. 



RESULTS 

FTS theory for QPTs. Just as the well-known 
regime of finite-size scaling in which the characteristic 
size L of a system is shorter than its correlation length 
£, an FTS regime arises when Td is the shortest time 
among the three time scales and thus dominates. The 
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theory of FTS [19J, |20| takes explicitly into account the 
rate R, which plays a role similar to L^ 1 since it imposes 
on a system an additional time scale that manipulates 
its evolution. For a length rescaling of factor b, an order 
parameter M transforms as 

M = b-^ v M{tb- z ,gb l ' v ', h z b^ lv , Tb z , Rb r , L~ l b), (1) 

where the two critical exponents ft and 5 are defined as 
usual in classical critical phenomena by M oc gP in the 
absence of an external probe field h z conjugate to M 
and M oc h 5 z at g = 0, respectively. Equation (JTJ in 
its classical form in which g is replaced by T — T c (T c 
is the classical critical point) can be derived from the 
renormalization-group theory (l9l . |20| . Here we apply it 
to the quantal case as critical singularities exhibit too 
albeit due to quantum rather than thermal fluctuations. 

With equation ([1} , one can describe in a unified frame- 
work different kinds of driven dynamics via changing g, 
h z or T and readily define different regimes and their 
crossovers. Taking for example g = Rt, neglecting h z 
and choosing b such that Rb r becomes a constant, one 
finds an FTS scaling form 



M = R p/ur h{gR~ 1/vr iTR- z ! r , L~ X R- 



l/r\ 



(2) 



where r 



19, 



X/v obtained from g = Rt and its rescaling 



20] and the function fa with an integer i denotes a 
scaling function. FTS dominates when \g\R~ 1 l vr <C 1, 
TRT z l r < 1 and L~ X R~ X I T < 1. The first is r d < 
\g\~~ vz ~ t s and the second Td <C 1/T = tt as they should 
be. Crossovers to other regimes occur near \g\ ~ R x / vr 
and T ~ R z l r as depicted in Fig. [T] and a similar one 
for L. The first gives \i\ ~ RT V z /^+ vz ) because g = Rt. 
This is just the result of the KZM 0, H in the 

thermodynamic limit (L — > oc) and at T = 0. 

Scaling for other observables can be simply written 
down. The density of quasiparticle excitations n ex , for 



instance, is 

n cx = R d l r h{gR^ lvr , TR- Z ' r , L^R- 1 ^), (3) 

as it is proportional to £~ d in a <i-dimensional space. Q 
gives rise to the leading singularity n ex oc R d l r ', which is 
the famous scaling of KZM or more precisely, 

n cx = R d / r f 2 {gR~ 1 / vr ) = R d/r f 2 (t/i) for T = in the 
thermodynamic limit in agreement with that suggested 
previously 0, It also predicts that n cx oc T d l z and 
ri cx oc L~ d with appropriate conditions in the QCR and 
finite-size scaling regime, respectively. Compa ring with 
the scaling forms in the imaginary-time quench [23J , here 
we directly focus on the real-time behavior which is real- 
izable directly in experiments. Different from similar ex- 
pressions including real-time dynamics [24], [25| and the 
finite temperature effects one sees that FTS pro- 
vides a unified framework not only to understand the 
nonequilibrium evolution, but also to derive various scal- 
ing forms. Indeed, we can readily consider other kinds of 
sweepings near the critical point as follows. 

Instead of sweeping g, when h z = R z t, one replaces R 
with R z and obtains similarly the order parameter 

M h = R%/ vr * x 
f3{gR- 1 '"',h,R-W"%TR-*/ r ',L- 1 R; 1 / r ') (4) 

with r z = z + ftS/v now. Different regimes and their 
crossovers can be readily defined in this case too. Dif- 
ferent from sweeping the quadratic coupling term, here 
we change the symmetry breaking field h z and we shall 
find this sweeping provides an experimentally realizable 
method to determine the critical point. 

We can also vary T with a rate Rt, i.e., T = Rrt, the 
scaling form for the order parameter is then 



j4(gn T ,n z H T ,1K T 



M T = R^ vrT x 



-z/r T j 

T 1 i I** 



-1 



K 



-l/r T \ 



(5) 



with tt — 2z. Theoretically, according to Feynman path 
integral representation varying T is equivalent to 
changing the imaginary time scale tt- Experimentally, 
varying temperature can even be realized in cold atoms 
experiments [26| . 

According to these scaling analyses, the regime dom- 
inated by Td, i.e., the FTS regimes can then be deter- 
mined. For changing g, a schematic phase diagram is 
shown in Fig. [TJ At T = 0, the FTS regime is thus 
just the impulse regime of the KZM. As a result, hys- 
teretic dynamics in this regime is expected to be de- 
scribable by FTS. One sees from Fig. [Q that the FTS 
regime pushes the QCR to higher temperatures, because 
only then tt dominates. Consequently, the FTS enables 
one to probe directly the quantum critical point and its 
scaling behavior at nonzero temperatures as the tem- 
perature becomes subordinate, i.e., just a perturbation. 
In other words, the nonequilibrium driving appears to 
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extend the zero-temperature critical point to the whole 
FTS regime, or conversely, has a 'dynamic cooling' ef- 
fect that enables one to probe zero-temperature scaling 
at nonzero-temperatures while keeping the temperature 
subsidiary. This offers us an extra approach to detect and 
study quantum criticality at finite temperatures. There- 
fore, FTS provides not only a unified, simple and feasi- 
ble approach to test predictions of microscopic theories 
but also a systematic approach to study equilibrium and 
nonequilibrium dynamics of quantum criticality. 

Determination of critical properties. Here we 
show that FTS can also provide methods to detect quan- 
tum critical properties such as the critical point and crit- 
ical exponents. As stated above, in the FTS region, ef- 
fects of the temperature can be neglected, so we consider 
T = in the thermadynamic limit L — oo. According to 
equation ([J}, the external field at Mh = 0, denoted by 
h zm , is scaled as 



h z0 (g,R z ) = R^higR- 1 '^). 



(6) 



At h z = 0, Mh is reduced to 

M (g,R z ) = Ri'^'higR- 1 '"'). (7) 

Differentiating M with respect to h z in (f4]), one obtains 
the susceptibility at zero field, 



X (g,R z ) = RM-Ww* f 7 (gR-y^y 
To fix the critical point, we can define a cumulant, 



C(g,R z 



h z oX 



fs(gR 



— 1/ VTz 



), 



(8) 



(9) 



similar to the Binder's cumulant in the finite-size scaling 
271 ] . As C is a function of only one independent vari- 
able, its curves of different R z intersect at the critical 
point h xc or g = at which C becomes a constant /s(0) 
independent on R z . 

After fixing the critical point, all the critical exponents 
can be estimated. For example, j58/vr z and j3/vr z can 
be estimated respectively from © and ([7]) by fitting h z0 
and M for a series of R z at g = 0. Similarly, from @ 
at T = and L = oo f3/vr can be estimated by fitting 
M for a series of R at g = 0. From these three exponent 
ratios and the scaling law [l[ [3(8 + 1) = (d + z)v, one can 
determine all the critical exponents. 

Application to the transverse-field Ising model. 
As an example of the FTS method to determine crit- 
ical properties, we consider the one-dimensional (ID) 
transverse-field Ising model whose Hamiltonian is 



A' 



71=1 



N-l 

£< 

n=l 



(10) 



and has been realized in CoNb206 experimentally [28], 
where cr^ and are the Pauli matrices, h x is the trans- 
verse field, and the Ising coupling has been set to unity as 



our energy unit. The model exhibits a continuous QPT 
from a ferromagnetic phase to a quantum paramagnetic 
phase at the critical point h xc (and so g = h x — h xc ) 
at T = [1]. The order parameter is the magnetiza- 
tion M — J2 n =i( a n) /N for the N spins with the angle 
brackets denoting the quantum and/or thermal average. 
As a method to probe the transition, we add to H a 
symmetry-breaking term —h z J~* * er^. 

We demonstrate our approach at T = at which some 
exact results are available for comparison. We solve 
the model using the time-evolving block-decimation al- 
gorithm [29j j . which is capable of treating large system 
sizes. We determinate the critical point in Fig. [2] and 
apply it purposely to determine the critical exponents in 
Fig. [31 The good agreement of the results collected in Ta- 
ble U shows the power of FTS. In Fig. [H these results are 
verified by scaling collapses, which provide also a method 
to estimate the range of the FTS regimes. Similar scal- 
ing collapses have been reported 13, 14, 23- 25|, however, 
here the critical exponents and the crossover regions are 
estimated by the FTS method self-consistently. 

Nonequilibrium FTS at finite temperations. 
Having successfully demonstrated FTS at T = 0, we now 
turn to nonzero temperatures at which most realistic ex- 
periments operate. In order to study thermal fluctuations 
and examine the general nonequilibrium FTS scaling Q , 
we have to consider an open many-body quantum system 
interacting with a heat bath [3(| . The state of system can 
be described by a density operator p according to quan- 
tum statistical physics. For weak system-environment 
couplings, after assuming Markovian and tracing over the 
bath variables, one obtains the master equation for the 
density matrix p in the Lindblad form [21 , , 



dp/dt — —i\H, p] + Cp, 



(11) 



where Cp = -c£S, iA ^{V^V^p + pV^V^ - 
^Vi^jpV^j)/2, c measures the coupling strength be- 
tween the system and the bath, Ne is the total number of 
the energy levels, ft = exp(- J B J /T)/Trexp(-W/T) with 
Ei being the «th eigenvalue of the H, and Vi-^j is the 
thermal jump matrix whose element at the jth row and 
ith column is one and zero otherwise in the energy rep- 
resentation. Vi^fj fulfills (3ipEVi-).j — PjVi^jPE, where 
Pe = cxp(— %/T)/Trexp(— H/T) is the equilibrium den- 
sity matrix operator whose eigenvalues are ft. This can 
be regarded as the detailed balance condition for Vi^j in 
equilibrium. 

The meaning of (jlip is quite clear: If the second term 
on the right hand side is neglected, it is the quantum 
Liouville equation contributing quantum fluctuations to 
the evolution of density operator p\ while if the first term 
on the right hand side is omitted, the diagonal part of 
the remaining equation is just a thermal master equa- 
tion contributing thermal fluctuations to p. The Lind- 
blad equation (jlip thus is a real-time dynamic equations 
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which integrates both the quantum and the thermal con- 
tributions. It has been widely used in quantum optics 



the rapidly developing numerical renormalization-group 
methods [37f , for example, seem quite promising to solve 



31] and relaxation processes in open quantum systems 

In order to show that nonequilibrium FTS does valid 
for T ^ 0, we solve numerically equation (1111) for the 
Hamiltonian (1101) along with the sweeping probe field by 
brute force. We find that the order parameter Mh can 
now saturate correctly with thermal fluctuations. More- 
over, the results displayed in Fig. [5] show clearly the va- 
lidity of the FTS scaling form (@}. 



DISCUSSION 

We have verified that the theory of FTS indeed pro- 
vides a unified description of the driving dynamics near 
the critical point from the realm of time. We have also 
determined the critical point and critical exponents accu- 
rately according to the scaling forms of FTS at T = 0. As 
mentioned, in the FTS regime, the effects of the temper- 
ature are suppressed by the nonequilibrium driving, i.e., 
thermal fluctuations or heating of the low-energy critical 
modes are governed by the driving time scale in the 
nonadiabatic nonequilibrium FTS regime. Only at higher 
temperatures do the thermal fluctuations take over and 
the QCR emerges. This 'dynamical cooling' effect allows 
one to determine the critical properties at finite temper- 
atures too. As experiments are always operated at finite 
temperatures, the FTS method can thus be used in the 
experiments directly. 

In QPTs, the Hamiltonian determines not only the par- 
tition function, but also the time evolution of any observ- 
ables. Further, the imagine time acts as z extra dimen- 
sions in the vicinity of a quantum critical point through 
the quantum-classical mapping. These place special im- 
portance on the dynamic critical exponent z. Note that 
z can be a noninteger for some frustrated systems [33j . 
Traditional methods employ usually a finite size in the 
imaginary-time direction and utilize finite-size scaling to 
estimate z 33]. In contrast, as the scaling forms shown, 
z enters the nonequilibrium scaling naturally. Thus FTS 
provides at least an alternative method to estimate z di- 
rectly from real-time scalings. 

We have used the Lindblad equation to describe the 
thermal behavior near the critical point. Note that it 
is a challenge to study the time evolution of nonequi- 
librium systems with many degrees of freedom even in 
isolated situations (3. IR l34l436| . Further, none of the an- 
alytic, semiclassical, or numerical methods of condensed- 
matter physics yields accurate results for dynamics in 
the QCR except for some special systems in ID Q. Al- 
though here we only solve directly the Lindblad equation 
(fTT]) for small lattices, the results show that it is suitable 
for describing the nonequilibrium behaviors at the finite 
temperature near the quantum critical point. Moreover, 
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the equation for larger lattice sizes 

In conclusion, FTS not only lights up the dark impulse 
regime of KZM at zero temperature but also sheds light 
on the QCR at nonzero temperatures by establishing its 
own regime. It also offers a powerful unified approach 
amenable to both numerics and experiments to study 
equilibrium and nonequilibrium dynamics of quantum 
criticality. We have also shown that the Lindblad master 
equation can be a valuable framework for studying real 
time quantum criticality at nonzero temperatures. Al- 
though we have studied a simple model for illustration, 
our approach should be applicable to more complex sys- 
tems as well. In addition, as our results indicate that the 
classical theory of FTS with appropriate modifications 
can well describe quantum criticality, new physics may 
be in action Q if it is violated. 



METHODS 

Solving the Schrodinger equation at zero tem- 
perature. The driving dynamics evolution at zero tem- 
perature is solved based on the time-evolving block- 
decimation algorithm [2!| . The efficiency of this method 
depends on the fact that singular values have an approx- 
imately exponential decay. After expressing the wave 
function in a matrix product state via Vidal's state de- 
composition, one can cut off the singular values with 
a prescribed error induced by the truncation. In the 
present work, we choose the lattice size of L = 2000 to 
eliminate the size effects. The number of singular values 
are kept at 16 which guarantees that the errors near the 
critical points are less than 10~ 5 . 

Solving the Lindblad equation at the finite tem- 
peratures. The Lindblad equation is solved by a finite 
difference method to second order with the time-interval 
St = 10~ 3 and the resultant errors are about 10~ 4 . 
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TABLE I: Critical point and exponents for the ID transverse- 
field Ising model 







P 


S 


v z 


Numerical 


0.999(2) 


0.125(11) 


14.9(6) 


0.98(4) 1.01(3) 


Exact [1] 


1 


0.125 


15 


1 1 




FIG. 1: Schematic phase diagram near a quantum crit- 
ical point under an external driving. Two equilibrium 
phases (green) in which the reaction time dominates are sepa- 
rated by two crossover domains fanning out from the quantum 
critical point O. One domain is the quantum critical regime 
QCR (blue), in which the imaginary time scale or temper- 
ature T dominates. It extends from the usual well-known 
QCR bounded roughly by T ~ ±\g\ uz at R = to finite driv- 
ing rates R. The other is the new finite-time scaling (FTS, 
red) regime in which the driving time dominates. At T = 0, 
it overlaps with the impulse regime enclosed roughly in the 
curves R ~ ±|<?|" r , thus improving the understanding of this 
dark regime of the Kibble-Zurek mechanism. The FTS regime 
extends to low T and pushes the QCR to higher T. Crossover 
occurs near T ~ R z ^ r . Because T is only subordinate in the 
FTS regime, the nonequilibrium FTS provides a valuable ac- 
cess at finite temperatures to the zero-temperature quantum 
critical point. 



0.98 0.99 1.00 1.01 1.02 



FIG. 2: Estimation of a quantum critical point. Curves 
of the cumulant C (Eq. ((9J) of different R z intersect at the 
critical point h xc or g = 0. Owing to possible errors from the 
truncation of the singular values in the Schmidt decomposi- 
tion [2S|, however, the intersections are slightly scattered as 
shown in the inset. Nevertheless, the average of all the inter- 
sections is h^c = 0.999(2), a good estimate of the exact value 
hxc = 1. We choose a lattice size of L — 2000, which is large 
enough for the size effect to be ignored. 




FIG. 3: Estimation of critical exponents. In the ther- 
modynamic limit and at T = 0, from the finite-time scaling 
forms 0, © and ©, M oc B?/ UT *, (a), h z0 oc R l3S/l ' r *, (b), 
and M oc R^l vr , (c), respectively, at the critical point g = 0. 
Our solutions with h* c = 0.999, T = and L = 2000 yields 
P/vr z = 0.0436, pS/ur z = 0.651 and fi/vr = 0.0622 from 
power-law fits. Using the scaling law !l] /3(S + 1) = (d + z)v, 
we then obtain all the critical exponents listed in Table[l]with 
their corresponding exact results for comparison. As the sta- 
tistical errors of the fits are tiny, we fit data at h xc — 0.997 
and 1.001 and the largest difference in each exponent is used 
as an estimate of the error given also in Table [I] 
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FIG. 4: Scaling collapses and crossovers. The original 
curves in insets (al) and (bl) collapse as expected when 
they are rescaled, because they are now functions of a single 
rescaled variable: Mh_RJ /3// " ri = f3(h z R z ~ l35 ^ rz ) at g = 0, 
(a), and MR~ l3/vr = h{gR- 1/,/r ), (b), according to the 
finite-time scaling (FTS) Q and respectively, at T = 
and large L. This therefore generalize the Kibble-Zurek scal- 
ing into the whole FTS regime. This regime, denoted by FTS 
and delineated roughly by the vertical lines in insets (a2) and 
(b2), crossovers to adiabatic evolutions. The exact positions 
are non-universal constants depending on the model. One 
sees that the FTS regimes are relatively large for the model. 




FIG. 5: Nonequilibrium scaling at nonzero tempera- 
tures. Data of Mo versus T plotted in the inset for the three 
different sets of R z and L so choosing as to fix the value 
of L~ 1 Rz 1 ^ T * collapse as expected onto a single curve de- 
scribed by M R- p/ur * = f 3 (TR7 z/rz ,L- 1 R7 1/rz ) for a fixed 
LR z ~ 1 ^ Tz = 0.777 according to the finite-time scaling Q at 
the critical point, c = 0.8 and h x = 0.999. 



